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ABSTRACT 

A new approximate theory which links the inherent flaw concept with the theory 
of crack tip stress singularities at a bi-material interface has been developed. 

Three assump tions were made: (1) the existence of inherent flaw (i.e., damage zone) 
at the tip of the crack, (2) a fracture of the filamentary ocmposites initiates at a 
crack lying in the matrix material at the interface of the matrix filament, (3) the 
laminate fails whenever the principal load-carrying laminae fails. This third 
assumption implies that for a laminate consisting of 0° plies, cracks into matrix 
porp^rv^imiar- to the 0° filaments are the triggering mechanism for the final failure. 

Based on t his theory, a parameter which is similar to the stress intensity 
factor for isotropic materials but wi th a different dimension was defined. Utilizing 
existing test data, it was found that Kq can be treated as a material constant. 

Based on this finding a fracture mechanics analysis methodology was developed. 

The analytical results are correlated well with test results. This new 
approximate theory can apply to both brittle and metal matrix composite laminates 
with crack or hole. 


INTRODUCTION 

The failure modes associated with fractures in fiber-reinforced ocmposites 
differ considerably from those of homogeneous isotropic materials. The failure modes 
in ccnposites are typically in the forms of transverse cracking, delamination, fiber 
breaks, matrix yielding, matrix cracking, fiber pull-out and fiber/matrix debonding. 
As a result, it is very difficult to predict the fracture behavior of composite 
materials exactly. 

In the past, a great deal of effort has been expended to investigate the 
fracture behavior of ocmposites. A number of fracture mechanics theories have been 
proposed. These theories have been reviewed and presented extensively in References 
1 and 2. 

* In this paper a microscopic theory which was originally proposed by Mar and 
Lin 9 has been modified into a new theory. This new theory links the inherent flaw 
concept^ , which postulates that a damage zone exists at the tip of crack, with the 
theory of crack singularities at a bi-material interface 9 . This combined theory 
can be to predict the notched strength of organic and metal matrix ccnposites 

with either a crack or a hole. 


♦This paper cites references 1-20. 
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THEORY 



The stress distribution at the crack tip in a thin plate for a homogeneous, 
isotropic elastic solid in terms of the coordinates shown in Figure 1 is qiven bv 
equation (1) 1 
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where 


on — gross nominal stress 

For an orientation directly ahead of the crack ( 0 = 0) 
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Irwin pointed out that equation (1) indicates that the local stresses near a 
crack depend on the product of the nominal stress and the square root of the half- 
flaw length. He called this relationship the stress intensity factor K, where for a 
sharp elastic crack in an infinitely wide plate, K is defined as 


K « On v/aJT 


(3) 


In the approach using linear-elastic fracture mechanics (TFFM) , k is a material 
parameter and may be determined from tests. 


For a finite-width plate, equation (3) is modified to 


K-Yon^ 


(4) 


Where Y is a parameter that depends on the plate and crack geometry. 

To develop similar concept for composite materials, the assumptions of 
references 3,9 were adopted in this paper; i.e. ; 
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. Hie existence of an inherent flaw (also called a damage zone) 
at the edge of a hole or at the tip of a crack. 

. Fracture of a filamentary composite initiates at a crack 
lying in the matrix material at the interface of the 
matriiQ/f ilament . 

. A laminate fails whenever the principal load-carrying laminae 
fails. This implies that cracks in the matrix perpendicular 
to the 0° filaments are the triggering mechanism for the 
final failure. 

Based on the above assumptions, the follcwing theories are developed: 

Using the same concept of stress intensity factor as is formulated above for 
isotropic materials, a material parameter similar to K is defined for composite 
material as 

"TT"= V on (ao) m (5) 

where m is the order of singularity of a crack whose tip is at the interface of two 
different materials as shown in Figure 2. Calculations for determining m are 
presented in Reference 18. 

Note that K has a different dimension from K. (K has a dimension of "stress 
tines length to the 1/2 power", while K has a dimension of "stress times length to 
the m power".) 

Although seme composite materials (such as polymeric matrix composites) fail in 
a brittle manner, a damage zone does develop which is analogous to the plastic zone 
for ductile materials. Using this concept in conjunction with equation (5) yields: 


—= Y o N (ao + C 0 ) m (6) 

where C Q is defined as an inherent flaw size. The term inherent flaw size is used 
since unnotched strength, oo , of a composite laminate is given by equation (6) for 
the case of vanishing 


K « °o (C 0 ) m 


( 7 ) 


where Y Q is the correction factor for infinite plate. 

It should be noted that C Q does not physically refer to an inherent crack, but 
a characteristic dimension of damage zone at the tip of a notch or crack prior to 
ultimate failure. 

The question we may ask new is whether K and C Q are material constants. Before 
we reach a conclusion, certain equations are helpful in answering this question. 
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Substituting equation (7) into equation (6) , and after seme manipulations, we 
obtain the following important equations that will be used to determine parameter""K7 
C 0 and notched strength of composite laminates 


a 0 


\ Op 

Vo N 


Vm 


- 1 


( 8 ) 


K 


air Y 0 Oo 




( 9 ) 


On _ 1 | Co | 

oo ^ [ao + Coj 

( 10 ) 

or 



( 11 ) 


In the following sections, K will be called the equivalent stress intensity 
factor for ccrposite materials. 

DENOMINATION OF EQUIVALENT STRESS INTENSITY FACTOR, IT 
AND INHERENT FLAW SIZE C Q 

Reference 8 provides extensive fracture test data of boron/aluminum 
laminates with various proportions of 0° and -45° plies. Hence, this test data will 
be used to characterize the fracture behavior of boroiy aluminum composite laminates. 
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FOtm/ATFhfr STRESS INTENSITY FACTOR. K 


From Reference 18, the order of stress singularity at the boron/alurainum 
interface, m equals to .347. Substituting m = .347 in equation (9), the equivalent 
stress intensity factor for boron/ aluminum composite laminate with center crack can 
be written as follows: 


K - ao 347 a 0 . 




( 12 ) 


Equation (12) is used to characterize the critical equivalent stress intensity 
factor of boron/aluminum laminates with various layups. 

By using the fracture test results from Reference 8, equivalent stress 
intensity factors were calculated from equation (12) and are tabulated on Tables 1 
through 4 for boron/aluminum composite with laminate constructions [0] 6 ,p [C>2/-45] s , 

(-45/0 2 ] s and [0/-45] respectively. Note the test results shown in Tables 1 to 4 
are average test results of Reference 8 . As shewn in Tables 1 to 4 , TT values seem to 
be a material property and vary with different laminate orientations . R values are 
also plotted on Figures 3 to 6. R AV G is the average value of K from sill the crack 
sizes. As shewn in the figures, can be approximately treated as a material 

constant. It has to be pointed out nere that "fcy G was obtained by averaging three 
different widths of plate. For w = 101 . eiratTK - values are almost the same for 
different ratios. 

Hie detailed calculations of TC are shewn in Reference 18 . 

Hie above tests are for center crack specimens. For other crack types and 
locations , the calculated equivalent stress intensity factors are shewn on Table 5, 
and are obtained from Reference 18 . It can be seen that "IT for these unidirectional 
boron/aluminum composites is constant for different crack conditions. Hie reason for 
the slightly different K as compared to Table 1 is due to a different value of 
ultimate tensile strength. 
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INHERENT FLAW SIZE. C Q 


Two methods were used to calculate the inherent flaw sizes for a composite 
laminate with center crack. 


Least Square Fit 

Equation (8) , in which C Q is a proportional constant, can be rearranged to yield 



By using the fracture test data as shewn in Tables 1 to 4, and the least square 
fit, C Q for various laminate constructions are determined as shown in Table 6. 

Average Equivalent Stress Intensity Method 

From Equation (7) we have 


Kavg = oo (Co) 


m 


(14) 


where m = .347 for borory'aluminuin. 


The inherent flaw size can be derived from equation (14) as follows: 
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In the case of boron/aluminum ccrposite, equation (15) becomes 


Co - 


Kavg \ 
°° , 


(16) 


The inherent flaw sizes for various laminate constructions were calculated using this 
method and were also tabulated on Table 6. 


We pointed out earlier that C Q does not physically refer to an inherent crack, 
but rather to a characteristic dimension of damage zone at a crack tip, prior to 
fracture failure. Comparing the dimension of C Q with respect to the crack size of 
B/Al specimens as shown in Table 1 to 4, it is clear that C Q cannot be neglected in 
the calculations of equivalent stress intensity factor. The significance of C Q will 
be further discussed later in this paper. 

Equation (7) is used to calculate critical equivalent stress intensity factor 
based on the C Q determined from least square fit method. These values are also 
plotted on Figures 3 to 6 as T^ qF . It can be seen from the plots that there is not 
much difference between am K™ except in the plot for the [0/+45] s laminates, 

where though there is a greater difference between ~K^y G and "^pspT^vG provides a 
better result. For this reason K AVG will be adopted in this paper and will be 
denoted as"!^. 


APPLICATIONS 

Once the critical equivalent stress intensity factor - ^ is known, the fracture 
strength of the composite laminates can be obtained from equations (10) and (11) . 

From equation (11) , we have a fracture strength prediction formula as follows: 


On 1 



1 + 


r Kq s 
k V.o° 



(17) 


From Equation (17) , it can be seen 
lesser the notch sensitivity and so from 
larger the inherent flaw size (i.e., the 
sensitivity. 


that the larger the term TVo 0 

the 

equation (10) , we can conclude that the 
damage zone size) , C Ql the lesser the notch 


In the following subsections, the theory developed here is used to predict the 
fracture strength of various ccrposite laminates. 
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BQFQi/AUMINUM (B/Al) OCT-fPOSITE 

Notched Strength Prediction 

For B/Al ccrrposite laminates, m = .347 and equation (17) becomes 
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(18) 

where for a composite laminate with a center crack, Y is assumed to be the same as 
that for an isotropic material 8 . 



For convenience, data from Tables 1 to 4 and 6 are summarized on Table 7 to be 
used for the following analysis. Substituting <j Q and T^, from Table 7 for different 
ply conditions into equation (18) values foroN/<j 0 can be obtained for various crack 
sizes and cure tabulated on Column 7 of Tables 1 to 4. The error shown on Column 8 is 
defined as the difference between calculated results and test results divided by test 
results. The analytical results are also platted on Figures 7 to 10. 

As can be seen, the prediction represents the experimental results reasonably 

well. 


ThefC values obtained from center cracked specimens are also used to calculate 
the strength of B/Al specimens with center holes. Comparison of the analytical 
results and test results 14 are shewn in Figures 11 to 13. Figure 11 shows excellent 
correlation between test and calculated results for [0]^ laminates with center 
holes, while in Figure 12 and 13, the maximum percentage error for analytical results 
is around 13%. The detailed calculations are in reference 18. This confirms the 
findings of reference 19 and 20 that the length of discontinuity and not the shape 
appeared to control the fracture strength. 

Notch Sensitivity of Boron/Aluminum (B/Al) Composite Taminat e 

The Ojj/oq test data in Tables 1 to 4 are plotted on Figures 14 to 16 for 
composite laminate of various widths to check the notched sensitivity of B/Al 
ccsrposites. It is obvious that [-45/0 2 ] s is the most notch-sensitive, while [0/-45] 
is the least sensitive to the notch size. For 2a c /W <.l [0] 6T is more sensitive to S 
the notch size than [0 2 /-45] s . For 2a c /W >.l, the C jj/oo vs ^a^/W curve for [0] 6T and 
[0 2 /-45 ] s laminates are almost identical. The trend mentioned above can be detected 
by using a parameter 15 defined by the ratioTu/ O0 as shown in Table 7. The ranking 
of various laminate configurations for theTCy/ O0 index is [0/-45] s , [0 2 /-45] , [0] 6T 
and [-45/0 2 ] s in descending order. The relative values of Xv/ CTq show a simifar trend 
to that of damage zone size, C Q . 
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It can be concluded that the larger the ratio Kq/oo (or the damage zone size 
C Q ) the less fracture sensitivity there is to the crack size. 

GRAPHITE /EPOXY fGr/Eth COMPOSITE 

Equation (9) is used to characterize the for Gr/Ep composites 6 . The 
detailed calculations are shewn in reference 18. TC for various ply orientations is 
plotted as shewn in Figures 17 to 19. It can be seen that Kq for Gr/Ep composites 
can be treated as a material constant. Table 8 summarizes tne characterization 
results of Gr/Ep composite laminates. For Thomel 300 graphite fibers in Narmco 5208 
epoxy resin, m = .297. 

Notched Strength Prediction and Notch Sensitivity 
For Gr/Ep composite with center crack, equation (17) becomes 


On 
o o 


Y 



- 3.367 


-.297 


3q 


(19) 

Substituting TJ1 and oo from Table 8, into equation (19), the fracture strengths 
of graphite/epoxy for various laminate constructions can be obtained and are plotted 
on Figure 20. The detailed calculations are shewn in reference 18. As can be seen 
from the figure, the correlation between calculated and experimental results is very 
good. 


It can also be seen from Figure 20 that the ratio TCy/ oo (or the inherent flaw 
size) in Table 8 can be used as a notch sensitivity indicator. The [0/90/-45] s 
laminate is less notch sensitive than the [0/-45] s and [0/-45] 2s laminates and 
accordingly it has a larger ratio of Kq/ 0q ( or c 0 ) than the other two laminates. 


COMPARISON OF ANALYTICAL RESULTS BETWEEN ANISOTROPIC AND NEW APPROXIMATE MODFTS 


Composite materials made by combining two materials with different elastic 
moduli are by nature anisotropic in the gross sense. The a nis otropic model for 
ccirposite materials is to assume that the oqrposite is a homogeneous, anisotropic 
solid. For an anisotropic fracture m = .5. Applying the inherent flaw concept 3 
for an anisotropic mode with center crack, we have 


Kq “ Yow^ao + Co) ^ 


( 20 ) 
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Kq - Oo(CS) 

(21) 

* 

Note that Kq has dimensions different from that of TC and C Q is the inherent 
flaw size corresponding to an anisotropic model. 

Examining equations (20) and (21) , we can derive the following useful equations 
for an anisotropic model 
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ao 



( 22 ) 
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(23) 


Kq “ oo 



( 24 ) 

Equation (22) can be used to obtain C* while equation (23) can be used to 
predict the fracture strength of composite laminates. 

Figures 21 to 24 shew the comparison of calculated results between the new 
approximate and anisotropic models. It is clear that the new approximate model 
predicts better results than the anisotropic model. 
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CONCLUSIONS AND RECCWMENDATICNS 


ODNCLDSIONS 


. Hie methodology developed here can be used to characterize the 
fracture toughness of the composite laminates and can be used as a 
design tool to predict the fracture strength of various composite 
laminates. 


. Hie parameter Kq which was called critical equivalent stress 
intensity factor is defined, and can be treated as a material 
constant for various composite laminates. 


. Hie new approximate model provides better results than those of the 
anisotropic model . 


. Hie larger the ratio — (or the damage zone size, C_., ) the higher 
the damage tolerance. °° 


RECOMMENEATIONS 


. Further verification of the new approximate theory with test results 
of various composite materials is needed. 


. Apply the theory developed here to predict the fracture strength of 
composite laminates with various crack angles. 


. Develop a methodology to predict the inherent flaw size at the crack 
tip before fracture. 
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Table 5. Critical Equivalent Stress Intensity Factor Of Unidirectional B/AI 5 Composite. 


SPECIMEN * 

2ao 2R 

or — 

w w 

Ravg MPa (mm) 347 

CH 

.25 

1196.00 

CS 

.40 

1227.00 

DEN 

.30 

1214.00 


* CH - CENTER HOLE SPECIMEN 
CS - CENTER SLIT SPECIMEN 
DEN - DOUBLE EDGE NOTCH SPECIMEN 


Table 6. Inherent Flaw Size, C 0 (mm). 



Least Square Fit 

Kavg Method 

[0]6T 

.633 

.552 

[02/^45] s 

.593 

.641 

[*45/02]s 

.405 

.457 

[0/±45]s 

.976 

1.28 





Table 7. Fracture Parameters For Various Laminate Configurations Of B/AJ. 


Ply 

Configuration 

0 0 

(MPa) 

Rq 

MPa (mm) 347 

Kq 

Oo 

MPa (mm) 347 

Co 

(mm) 

[0]6T 

1672 

1360 

.81 

.552 

[02/±45] s 

800.1 

685.6 

.867 

.641 

[±45/02] s 

910.5 

693.8 

.762 

.457 

[0/±45]s 

581.4 

633.5 

1.09 

1.28 


Table 8. Fracture Parameters For Various Laminate Configurations Of Gr/Ep. 


Ply 

Configuration 

Oo 

(MPa) 

Kq 

MPa (mm)' 297 

Kq 

Oo 

(mm)* 297 

[0/±45]2s 

541.0 

408.6 

.755 

(0/±45] s 

541.0 

393.5 

.727 

[0/90/±45] s 

454.0 

437.7 

.964 
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Figure 3. Equivalent Stress Intensity Factor For [OjeT B/AI Composite. 
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Figure 4. Equivalent Stress Intensity Factor For [02At45]s B/AI Composite. 
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Figure 5. Equivalent Stress Intensity Factor For [±45/02] s B/AI Composite. 
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Figure 6. Equivalent Stress Intensity Factor For [0/±45] s B/AI Composite. 














Figure 12 Notched Strength Prediction for B/Ai [02/*45]s Composites With Center Hole 
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Figure 13 Niched Strength Prediction For B/A! |0^4i»] s Composite With Cc-ntei Hole. 
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Figure 14 Notch Sensitivilies For Various B/A! Laminate Orientations (w = 19 1mm). 
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Figure 21. Comparison Of Analytical Results Between Anisotropic And New Approximate Models B/Al {0)^1 



Figure 22. Comparison Of Analytical Results Between Anisotropic And New Approximate Models B/Al (<>2/±45) s 






Figure 23. Comporison Of Analytical Result* Between Anisotropic And New Approximate Models B/AI (±45/0 2 ) 



Figure 24. Comparison Of Analytical Results Between Anisotropic And New Approximate Models B/AI (0/±45> s 





